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We present a measurement of the mass difference between top {t) and anti-top it) quarks using 
tt candidate events reconstructed in the final state with one lepton and multiple jets. We use 
the full data set of Tevatron -/s — 1.96 TeV proton-antiproton collisions recorded by the CDF II 
detector, corresponding to an integrated luminosity of 8.7 fb~^. We estimate event-by-event the 
mass difference to construct templates for top-quark signal events and background events. The 
resulting mass difference distribution of data compared to signal and background templates using 
a likelihood fit yields AMtop = Mt - Mj = -1.95 ± 1.11 (stat) ± 0.59 (syst) GeV/c^ and is in 
agreement with the standard model prediction of no mass difference. 

PACS numbers: 14.65.Ha, 13.85.Ni, 13.85.Qk, 12.15.Ff 



The laws of the standard model of particle 
physics (SM) are invariant under the simultaneous trans- 
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formations of charge conjugation, parity, and time rever- 
sal (CPT). Conservation of CP T is, therefore, fmidamen- 
tal and provides one of the most important constraints 
on the SM. However, examining any possibility of CPT 
violation is important, as there are well-motivated ex- 
tensions of the SM allowing for CPT symmetry break- 
ing li;]. In CPT- conserving models, particles and their 
anti-particles must have identical masses and widths. 
Thus, any mass difference between a particle and its anti- 
particle would indicate a violation of CPT. CPT invari- 
ance has been tested for many elementary particles such 
as leptons and hadrons [2. ^] , but not in the quark sector 
except for the top quark Q . For all quarks except the top 
quark, direct mass measurements are nearly impossible 
because the quark hadronization time scale is approxi- 
mately an order of magnitude less than the quark de- 
cay time. After hadronization occurs, only the masses of 
hadrons are observable and give, at best, only an approx- 
imate estimate of the constituent quark's masses. On the 
other hand, as the lifetime of the top quark is of the or- 
der of lO"^"* seconds, it decays before hadronizing and a 
precision measurement of its mass and of the difference 
between the quark and anti-quark masses can be made. 

Since the top quark discovery, close to three thousand 
of tt events has been collected at the Tevatron pp collider. 
This sample makes measuring the top-quark mass (Mtop) 
possible to an accuracy of approximately 0.5% (Mtop = 
173.2±0.9 GeV/c2) [sj and the mass difference (AMtop = 



4 



Mt — Mf) between t and t quarks to a comparable pre- 
cision. The DO collaboration performed a few measure- 
ments of AMtop using matrix element analyses 01 ■ 
The most recent DO result, based on a 3.6 fb~^ data sam- 
ple, reports AMtop — 0.8 ± 1.9 GeV/c^, consistent with 
zero as predicted in the SM. The CDF collaboration per- 
formed a measurement using a 5.6 fb~^ data sample [sj 
and found AMtop = — 3.3±1.7 GeV/c^ which is also con- 
sistent with zero to within two standard deviations. To 
date, the most precise measurement is performed by the 
CMS collaboration, AMtop = -0.44 ± 0.53 GeV/c^ @. 

This paper reports on the final CDF measurement of 
AMtop based on the full Run II data set corresponding to 
an integrated luminosity of 8.7 fb~ . We reconstruct the 
mass difference between t and t quarks in each data event 
and compare its distribution with template distributions 
derived from Monte Carlo (MC) model simulations to 
estimate AAitop- This is an update of a previous mea- 
surement that used a subset of the present data ^SJ. In 
addition to the larger data sample, we improve the jet 
energy calibration by applying an artificial neural net- 
work to achieve better jet energy resolution [10], as in a 
recent measurement of Mtop 111]. We also re-examine the 
systematic uncertainties resulting from the larger size of 
the control samples. 

In the SM, t and i quarks decay almost exclusively 
into a W boson and a bottom quark {t bW~^ and 
t — ^ bW~) [l3|- The case where one W boson decays 
to a charged lepton (electron or muon) and a neutrino 
(W+ £~^h' or W~ £~h' including the cascade decay 
of W —i' Ti' and r — > iv) and the other to a pair of jets, 
defines the lepton+jcts channel. To select tt candidate 
events in this channel, we require one electron (muon) 
with Et > 20 GeV (pr > 20 GeV/c) and pseudorapidity 
1 77 1 < 1. 1 I l3l| . We also require large missing transverse 
energy [Hi {^t > 20 GeV) and at least four jets. Jets 
are reconstructed applying a cone algorithm with radius 
AR = ^(Ar/)2 -I- (A0)2 = 0.4 [ll_. Besides the stan- 
dard jet energy scale corrections |16j . we use an artifi- 
cial neural network that includes additional information 
to the calorimeter one, such as jet momentum from the 
charged particles inside the jet [13] ■ This additional in- 
formation improves the resolution on the reconstructed 
jet variables, resulting in approximately a 10% improve- 
ment in statistical precision. Jets originating from b 
quarks are identified (tagged) using a secondary vertex 
tagging algorithm (17j . In order to optimize the back- 
ground reduction and to improve the statistical power of 
the measurement, we divide the sample of ti candidates 
into subsamples with zero (0-tag), one (1-tag), and two 
or more (2-tag) fe-tagged jets. 

For the 0-tag events, we require exactly four tight 
jets (transverse energy Et > 20 GeV and I77I < 2.0). In 
case of the 1-tag and 2-tag events, three tight jets and one 
or more loose jets {Et > 12 GeV and I77I < 2.4) are re- 
quired. To reduce background contributions in the 0-tag 
or 1-tag events, we require the scalar sum of transverse 
energies in the event, Ht = £^T^*°" + ^T + I]tour icts ^t*' 



to exceed 250 GeV. The Ht requirement is not applied 
to the 2-tag events because of the small background con- 
tribution in this subsample. We divide the 1-tag and 
2-tag samples into subsamples based on the number of 
tight jets. We denote as tight subsample the sample re- 
quiring exactly four tight jets and loose subsample the 
sample consisting of the remaining events. This results 
in five subsamples: 0-tag, 1-tagL, 1-tagT, 2-tagL, and 2- 
tagT, where T and L denote tight and loose subsamples, 
respectively. 

The primary sources of background events are VF-|-jets 
and QCD multijet processes. To estimate the contri- 
bution of each process, we use a combination of data- 
and MC-based techniques described in Rcf. [18]. For the 
Z-|-jets, diboson, single top quark, and tt events we nor- 
malize the number of simulated events using their the- 
oretical cross sections 19 21j]. We use the data-driven 
techniques described in Ref. [2^ to estimate the QCD 
multijet background. The T^-l-jets background shape is 
modeled using MC generated samples but the number of 
events is derived from the data sample by subtracting all 
other contributions, including the tt signal, from the data 
events. Table U summarizes the data sample composition. 

We assume that all selected events are lepton-t-jets 
tt events and reconstruct AAftop, event-by-event, using 
a special-purpose kinematic fitter Measured four- 
vectors of the lepton and jets are corrected for known ef- 
fects as described in Ref. fl&\ , and appropriate resolutions 
are assigned. The unclustered transverse energy (Ut) is 
estimated as a sum of all transverse energy in the calo- 
rimeters that is not associated with the primary lepton 
or with one of the leading four jets. It is used to cal- 
culate the neutrino transverse momentum. The longi- 
tudinal momentum of the neutrino is a free parameter 
which is effectively determined by the constraint on the 
invariant mass of the leptonically-decaying W boson. To 
estimate AMtop, we define a kinematic function, 

X = Si=£4jots(PT -Pt ) /^t 

I /TT^^ TTmca.s\2 / ^2 

+ ^k=x,v[UT^ - Ut^ ) /(Tfe 

+ (M,, - Mw f/Tl^ + (M,, ~ Mw)VT^w 

+ {M,j, - (M-- + dm,eco/2)}Vr? 

+ {Mu. - (ACp- - dm„co/2)}Vr?, (1) 

where dmreco is obtained at the lowest ^^id rep- 
resents the reconstructed mass difference between the 
hadronically- and leptonically-decaying top quarks, 
Mhjj — Mmv In Eq. (1), we constrain the lepton pT 
and the four leading jets pt to their measured values 
and uncertainties {oi). We also constrain Ut in the sec- 
ond term of Eq. [1] In the remaining terms, we constrain 
the W boson mass {Mw) to Mh^=80.4 GeV/c^ and the 
average of t and t masses to Mf™-=172.5 GeV/c^. The 
quantities Mjj , Mgi, , Mjjjj , and Mmh refer to the invari- 
ant masses of the particles denoted in the subscripts. The 
total widths of the W boson, Y\y = 2.1 GeV, and of the 
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TABLE I: Expected and observed numbers of sij 


;nal and backf 


ground events assuming 


a tt production cross 


section att — 


7.45 pb and i\ftop = 


172.5 GeV/c . 












0-tag 


1-tagL 


1-tagT 


2-tagL 


2-tagT 


W+jets 


778±219 


197±69 


114±42 


11.4±4.9 


8.0±3.4 


Z+jets 


55.7±4.9 


10.3±1.2 


6.7±0.8 


0.8±0.2 


0.5±0.1 


Single top 


5.1±0.4 


11.7±1.0 


7.2±0.6 


2.2±0.2 


1.7±0.2 


Ui boson 


Dd.9±5.9 


11.7±1.5 


9.0±1.2 


0.9±0.2 


0.9±0.2 


QCD niultijet 


133±107 


31.7±1.2 


20.9±16.9 


4.3±4.3 


2.9±3.5 


Total background 


1038±244 


262±70 


158±45 


19.5±6.5 


14.0±5.0 


tt signal 


620±83 


694±87 


847±105 


188±29 


294±45 


Expected 


1658±257 


957±111 


1005±114 


208±30 


308±45 


Observed 


1712 


919 


1018 


214 


286 



top quark, Tf — 1.5 GeV, are taken from Ref. [l3|. We 
assume that the total widths of the t and t quarks are 
equal. Determining the reconstructed mass difference of 
t and t, Am™'^", requires the identification of the particle 
type (t or i), which is achieved using the electric charge 
of the lepton (Qiopton), Amf = -Qiopton ■ drrircco- In 
the events with a positive (negative) lepton, t (t) decays 
leptonically and t (t) decays hadronically. Because of the 
different resolutions of the jets, lepton, and unclustered 
energy, the Amj°'^° distribution from the positive lepton 
events is different from that of the negative lepton events. 
To improve the resolution of the Aml'"^°, and allow using 
the appropriate distribution in the hadronic-to-leptonic 
and in the leptonic-to-hadronic mass difference, we di- 
vide each subsample into the two new subsamples based 
on the lepton charge. We then have ten subsamples in 
total. 

Assuming that the leading four jets in any event come 
from the four final quarks at the hard scattering level, 
there are 12, 6, and 2 possible jet-to-quark assignments 
for 0-tag, 1-tag, and 2-tag samples, respectively. The 
minimization is performed for each jet-to-quark as- 
signment, and AmY"^° is taken from the assignment that 
yields the lowest (Xmin)- The 6-tagged (zero 6-tag) 
events with xLin > 9.0 (Xmin > 3.0) are rejected due 
to the poorly reconstructed kinematical properties. To 
increase the statistical power of the measurement, we 
employ an additional observable, AmJ'^'^"'-^-', which cor- 
responds to the 2nd lowest in the jet-to-quark combi- 
natorics. Although it has a poorer sensitivity, AtoJ°'^°*-^'' 
provides additional information on AMtop and reduces 
the statistical uncertainty by approximately 10%. We 
use two observables {AmY"^° and Atoj°™^^'') simultane- 
ously for the measurement. 

Using MADGRAPH [^^l, we generate tt signal samples 
with A A/top between -20 GeV/c^ and 20 GeV/c^ in 
2 GeV/c^ intervals. Parton showering of the signal events 
is simulated with pythia [25], and the CDF II detector 
is simulated using a GEANT-based software package [26| . 

We estimate the probability density functions (PDFs) 
of signal and background using the kernel density es- 



TABLE II: Summary of systematic uncertainties on AAitop. 



Source 


Uncertainty (GeV/c^) 


Signal modeling 


0.14 


Parton showering 


0.17 


b and b jets asymmetry 


0.38 


Higher-Order effect 


0.16 


Jet energy scale 


0.07 


Parton distribution functions 


0.12 


6-jet energy scale 


0.05 


Background shape 


0.20 


Gluon fusion fraction 


0.05 


Initial and final state radiation 


0.10 


Finite Monte Carlo samples 


0.07 


Lepton energy scale 


0.06 


Multiple hadron interaction 


0.05 


Color reconnection 


0.23 


Total systematic uncertainty 


0.59 



timation [l^, [1^. We construct the two dimensional 
PDFs that account for the correlation between Am™'^° 

and Am].'^'^°^^\ First, at discrete values of AMtop from 
-20 GeV/c2 to 20 GeV/c^, we estimate the PDFs for 
the observables from the above-mentioned madgraph tt 
samples. We interpolate the MC distributions to find 
PDFs for arbitrary values of AMtop using the local poly- 
nomial smoothing method f29| . Then, we fit the signal 
and background PDFs to the unbinned distributions ob- 
served in the data using a maximum likelihood fit [30|. 
Separate likelihoods are built for the ten subsamples, and 
the overall likelihood is obtained by multiplying them to- 
gether. 

We test the fitting procedure generating 3000 sim- 
ulated experiments for each of eleven equally-spaced 
AMtop values ranging from —10 GeV/c^ to 10 GeV/c^. 
The fit estimates and their uncertainties in the simulated 
experiments are found to be unbiased. 

We examine a variety of systematic effects that could 
affect the AMtop measurement. To estimate the system- 
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atic uncertainties, we compare the results from simulated 
experiments in which we vary relevant parameters within 
one standard deviation. All systematic uncertainties are 
summarized in Tabic |TT1 The dominant source of system- 
atic uncertainty is attributed to a possible difference in 
the detector response between b and b jets. To estimate 
this effect, we select a bb sample by requiring exactly 
two 6-tagged jets per event using a sample triggered on 
jet {Et > 20 GeV) . In addition, one 6-tagged jet is re- 
quired to contain a soft muon from leptonic decay so that 
the charge of the b quark associated with the jet can be 
determined. The energy scale of b and b jets in the data 
is compared with di-jet MC events to estimate indepen- 
dently the difference of the b and b jet energy scale. The 
result (—0.44 ± 0.40%) is consistent between data and 
MC events. We perform simulated experiments by vary- 
ing the b and b energy separately. We estimate the sig- 
nal modeling uncertainty by using simulated experiments 
with events generated with madgraph and pythia. We 
also estimate a parton showering uncertainty by a pply ing 
different showering models (pythia and HERWIG [3l| ) to 
a sample generated with ALPGEN [s^I- Higher-order ef- 
fects are estimated using MC@nlo generator [s^. The 
background shape systematic uncertainty accounts for 
the variation of the background composition. The color 
reconnection systematic uncertainty [34i is evaluated us- 
ing samples simulated with and without color reconnec- 
tion effects in pythia tunes |35j . We vary the parameters 
of parton distribution functions to account for systematic 
effects. The jet energy scale uncertainty, the dominant 
uncertainty in most of the Mtop measurements, is par- 
tially canceled in the t and i mass difference. Other 
sources of systematic effects, including uncertainties in 
gluon radiation, multiple hadron interaction, finite size 
of MC samples, 6-jet energy scale, and lepton energy 
scale, give small contributions. The total systematic un- 
certainty of 0.59 GeV/c^ is calculated as a quadrature 



sum of the listed uncertainties. 
The resulting mass difference is 

AA/top = -1.95 ± 1. II (stat) ± 0.59 (syst) GeV/c^ (2) 

Figure [1] shows the observed distributions of the observ- 
ables used for the AMtop measurement. The density es- 
timates for tt signal events with AAftop = GeV/c^ and 
for background events are overlaid. 

In conclusion, we examine the mass difference between 
t and i quarks in the lepton+jets channel using CDF 
II data corresponding to an integrated luminosity of 
8.7 fb^^ bompp collisions at ^/s — 1.96 TeV. We measure 
the mass difference to be AMtop = Alt — = — 1.95 ± 
l.II (stat) ±0.59 (syst) GeV/c^ = -1.95 ±1.26 GeV/c^. 
This result is consistent with A A/top = GeV/c^ and 
conservation of CPT symmetry. 
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